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Abstract 

An approximation is used that permits one to explicitly solve the 
two-point Schwinger- Dyson equations of the U{N) lattice chiral mod- 
els. The approximate solution correctly predicts a phase transition for 
dimensions d greater than two. For d < 2, the system is in a single dis- 
ordered phase with a mass gap. The method reproduces known N = oo 
results well for d = 1. For d = 2, there is a moderate difference with 
N = oo results only in the intermediate coupling constant region. 
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I. Introduction 

The generation of a mass gap in two-dimensional spin systems is believed to be 
analogous to the generation of a non-zero string tension in four-dimensional non- 
abelian gauge theories [jl|, 0. Of particular interest are the d = 2 matrix chiral 
models. They are asymptotically free and have properties similar to d = 4 gauge 
theories [|I], ^ A solution to the N = oo chiral models has not yet been found. In 
contrast, it is straightforward to solve the 0{N) vector spin models as N —>■ oo. 

The degrees of freedom of the 0{N) models are iV-component vectors v{x) as- 
signed to the sites x of a lattice satisfying the constraint f (x) -f (x) = 1. If it were not 
for this constraint, the 0{N) vector model would be a theory of non- interacting 
particles. Thus the effect of the constraint is to introduce interactions. Due to these 
inter-particle forces, an interesting picture of mass generation arises [Q]. In terms of 
bare quanta, the interactions are repulsive and of order The quantum vacuum 

involves fluctuations in which any one of the quanta are created and destroyed. 
From a Euclidean viewpoint, such events correspond to particle loops, and so the vac- 
uum is like a gas of closed loops. Such loops are small and dilute in strong couping 
since the mass of a particle is large. As one moves to weaker couplings, the mass 
decreases and the loops become larger and more plentiful. If interactions could be 
turned off, the mass would eventually vanish and a phase transition would occur. 
For dimensions < 2, two large closed loops generically intersect. Because there are 
repulsive forces at such intersection points, a particular loop feels a pressure from the 
gas of surrounding loops. For d <2 and A^ > 3, this pressure prevents any particular 
loop from becoming too large and generates an interaction-induced mass which never 
vanishes. By these means, the vector models avoid a phase transition when d < 2. 
The models remain in the strong coupling phase in which the particles are massive 
and the system is disordered: Correlation functions fall off exponentially with dis- 
tance. In the large A^ limit, the vector model is exactly solvable: after interactions 
are incorporated into an effective coupling, particles are free. For d < 2, the effective 
coupling never attains a value sufficient to produce a phase transition to a massless 
spontaneously broken phase. 

The picture for the matrix chiral models is qualitatively the same [Q. Forces 
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at the intersection of particle loops are repulsive, as can be seen from Figure 6 of 
ref.0. However, no effective free-theory formulation of the model has been derived, 
although an exact 5- matrix has been proposed For SU{N), the spectrum is 

Mn = Msin (mr/N) /sin {tt/N), for = 1, 2, . . . , — 1. Interestingly, as ^ oo, 
the S'-matrix becomes the identity matrix to leading order in a expansion. This 
suggests that perhaps the large limit of the SU{N) chiral model is a free theory of 
particles. 

II. The Approximate Large N Method 

In this work, we obtain an approximate solution for the chiral models. Correlation 
functions and the mass gap are obtained. Because the method is similar to the one 
used in the vector models, it is useful to quickly review the vector models in the large 
N limit. The action A for the d-dimensional U{N) vector model is 

A = (3Nj2v*{x)-v{x + A) , (1) 

x,A 

where v{x) is an A^-dimensional complex vector satisfying v*{x) ■ v{x) = 1, (3 is 
proportional to the inverse coupling constant, and the sum over A is over the 2d 
nearest neighbor sites, i.e., A is ±ei, ±62, . . ., ie^, where is a unit vector in the 
ith direction. The U{N) model is equivalent to the 0{2N) vector model. 
The lattice Schwinger- Dyson equation for the two-point function is 

6.y = {v* (x) ■ V (y)) (^) ■viy + A)) + 

A 

(3j2{v*i^)-Hy) v*{y + A)-v{y)) . (2) 

A 

Since the four-point function enters in the last term, the vector model is not a free 
theory. However, large factorization allows the four-point function to be expressed 
in terms of two-point functions via 

{v* {x) ■ V {y) v*{y + A)-v (y)) = 

{v*{x)-viy)){v^y + A).v{y))+0^^'^ . (3) 
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The nearest neighbor expectation, {v* (y + A) ■ v (y)), is independent of y. A hnear 
equation for {v* (x) ■ v (y)) is thus obtained when is substituted in (^. Conse- 
quently, the N = oo vector model is a free theory governed by an effective coupling 

^'^^ " l + 2dp {v* {0)-v{A)) ■ 

For d < 2, no phase transition is encountered as /5 oo: The mass gap is a decreasing 
function of f3eff but the effective coupling does not reach the critical value at which 
the mass gap vanishes. The theory remains in the strong-coupling disordered phase. 
For d > 2, a phase transition occurs to a weak-coupling spin- wave phase in which 
Goldstone bosons appear due to the breaking of the global U{N) symmetry. This is 
in agreement with expectations [|, |^ . 

The action A of the matrix chiral model is 

A = pNj2Tr(u^x)U{x + A)) , (5) 

where U (x) is a unitary matrix in a group G. The chiral model is the non-linear 
effective low-energy theory of the massless quark model involving flavors. For 
this reason when d = 4, G = SU{N) is often used with A^ = 2 or A^ = 3. The 
matrices U{x) can be regarded as the Goldstone bosons of spontaneous broken chiral 
SUl{N) X SUr{N). 

In what follows, we treat the case G = U{N). The Schwinger-Dyson equation for 
the two-point function is 



S^y = 1 (Tr (^t (^) u (y)) )-f3^1.(Tr (t/t {x) U {y + A))) + 

f^T.lf{Tr (U^ (x) U {y) {y + A)U {y)) ) . (6) 

A 

Equation (||) is not useful in weak coupling where (3 becomes large because the last 
two terms individually become sizeable. Since the left-hand- side remains constant, a 
delicate cancellation between the two takes place. 

Equation (@) is similar to eq.(|D, except that matrix products arise instead of 
dot products. Consequently, the four-point function which enters in (^) cannot be 
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factorized in the large N limit. However, one may consider the following vector- 
model-like replacement 

1 (Tr (f/t (x) U (y) {y + A)U (y)) ) ^ 

aj^(Tr{uHx)U{y + A))) + bj^(Tr{uHx)U{y))) , (7) 

where a and b are arbitrary functions of (3, whose choices are at one's disposal. 

It is important to use values of a and b which reproduce the leading orders of 
perturbation theory, since one is interested in taking the continuum limit for which 
the coupling g goes to zero. Writing U {y) = exp [z$ {y)] and expanding about the 
identity matrix, one finds U (y) W {y + A)U (y) = 1 + 2z$ (y) - z$ (y + A) + . . .. 
Hence, if U (y) W {y + A)U (y) ^ all {y + A) + bU (y) then one needs. 



a = -l + 0[g') , b = 2 + 0[g' 
Alternatively, one can perform an operator product expansion. Write 

U (y) {y + A)U (y) = 2U (y) - f/ (y + A) + 

{U{y + A)-U (y)) W {y + A) (U (y + A) - U (y)) . (9) 

Note that the third term is a higher dimensional operator since it involves two deriva- 
tives. One again concludes that a and b have the perturbative expansions of eq.(|). 
After the substitution in eq.(|^) is made, one obtains 

1 {Tr (C/t (X) U (y))) = j^^^^ + A„ ^ 1 (Tr {U> (x) f/ (, + A))) , (10) 



P.tf= , . ■ (11) 



where 

(3{l-a) 
1 + 2dbf3 

One advantage of the substitution in eq.([7p is that the effective coupling /?e// remains 
finite as /5 — > oo so that the corresponding equation for the two-point function does 
not involve the above-mentioned delicate cancellation between large terms. 
One possibility is to choose 

a = -1 , 6 = 2Ga , (12) 
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where the average hnk is 

GA = ^(Tr(f/t(y + A)f/(y))) . (13) 

The choice in eq.(|T2|) produces the correct result in eq.(|^ when y = x. As a side 
remark, selecting a = and b = Ga yields the two-point function of the vector model 
in the large N limit. 

The the equation for two-point function G{x,y) = (Tt (u'' (x) U {y)j'^ / N in 
eq. ([ToD is linear. The solution is 

^ J exp[ip-{y-x)] 
-TT V J I- 2l3eff E cos(pi) 

i=l 

If one requires that G (x, y) = 1 ioi y = x, then the self-consistency condition 

" ' (15) 



l + 2df3bJ (27V)\ ^ f 



eff 

1=1 

must be satisfied. 

A more detailed analysis reveals that the substitution in eq.(^ with a and b 
given in eq.(|T^) is guaranteed to reproduce the 0(1) and O (g^) terms of correlation 
functions when a perturbative expansion of the theory is performed. Furthermore, 
eqs.(0) and (|T3|) lead to qualitatively the same results as in the vector model case: 
For d < 2, the chiral model is in the disordered phase and there is a mass gap for all 
values of j3. For d > 2, a transition occurs at a finite value of /3 to a phase in which 
global symmetries are spontaneous broken. This is in agreement with expectations 

i Si- 

III. The One-Dimensional Case 

When d = 1, the chiral model can be solved exactly. Solutions are known for the 
continuum and lattice cases at finite and for N = oo @, §]. This case allows one 
to test the approximate large A^ method. When d = 1, solving eqs.(^ and ([T^ gives 
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where 



A = ^^^=/3e// . (17) 

1 + 2/56 ^ ' 



As a consequence of the self-consistency condition in eq.(|T5D, one also has 



2(3 ^ ' 

For large (3, the mass gap m is 

mA~^ = — — ^ + . . . , (19) 

2{l-a)p 48/?2(l-«) 

where is the lattice spacing, which is often denoted by a. 

The continuum limit is obtained by taking /? ^ oo and A^^ with f3A = 
1/ (g^N) fixed. Here, gc is the continuum coupling constant. Hence, 

m = ^ . (20) 
4 

This is the exact d = 1 result for the U{N) gauge theory. 

IV. The Two-Dimensional Case 

When d = 2, the approximate large N method leads to a mass gap which is 
exponentially small as /5 — > cxd: 

mA-^ ^exp[-2TTp{l-a)] . (21) 

For a = —1, m ^ A exp [— 47r/3]. It is remarkable that an exponentially suppressed 
mass gap is obtained, given that the substitution in eq.(0) is guaranteed only to 
reproduce perturbative results to order g"^. However, the mass gap in eq.(|2l|) for 
a ^ — 1 is not as suppressed as the result obtained from perturbative renormalization 
group analysis, which gives m ^ Aexp [— Svt/S] | |10| |. 

One possibility to overcome this discrepancy is to permit the coefficients a and b 
in eq.(^ to depend on the distance between x and y. For weak coupling, a ~ — 1 at 
short distances but at large distances a ~ —3. However, one would have to find some 
theoretical justification for this dependence of a on \y — x\. 

V. The Matching of Both Strong and Weak Couphng 
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Although the values of a and b in eq. (|T2D produce good results for weak couplings, 
they fail to do so in the strong coupling region: For example, the mass gap behaves 
as —ln{2(3) + 0{(3) instead of —ln{j3) + 0(/3). To obtain results which are accurate 
in both strong and weak coupling, we have found that 

a = - (Ga)' , b = GA + (Ga)' (22) 

works reasonably well. These values of a and b satisfy the weak coupling limits in 
eq.(||) since Ga = ^—g^N/ {4:d)+0{g'^). In addition, one can show that they reproduce 
correlation functions correctly to O {g'^) in a perturbative expansion. It turns out that 
eq.(p2D produces the correct N = oo strong coupling expansions for the mass gap and 
Ga to order P^. Since eqs.(0) and (p^) give good results in both the strong and 
weak coupling limits, there is a reasonable chance that results are good throughout 
the entire coupling constant region. The idea of weak-strong coupling interpolation 
recently was used successfully in an approximate analytic computation of the 
glueball mass in the SU{3) gauge theory |[TT |. 

To examine how well eq. (|22|) works for intermediate couplings, one can used the 
solvable d = 1 C3iSG clS db guide. When N = oo, Ga = P ioi p < 0.5 and Ga = 
1 — 1/(4/3) for P > 0.5 H, |], and the two-point function is given by 

G(x,i/) = (Ga)I^-^-| . (23) 

Using the approximate large N method with eqs.(|T8|) and (^), one obtains the 
following equation for Ga: 

(Ga)^ + ^ = 1 ■ (24) 
For the approximate large N method, the two-point correlation function is again 
given by eq. (^3|) but with Ga determined from eq . (^4|) . Hence the difference between 
the exact and approximate results is governed by Ga. Figure 1 displays Ga. The 
approximate large N method yields values of Ga which are slightly less than the 
exact N = oo value. The difference reaches a maximum of 6.7% near P = 0.6. Very 
good agreement is, of course, obtained for P small or large. 

The lattice N = oo d = 2 chiral model has not been exactly solved, but Monte 
Carlo data is available for finite values of N and there exist analytic weak and strong 
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coupling results |12, 13) In Figure 2, we compare the approximate large N method 
for the values of a and b in eq.(^) versus weak and strong coupling results. For the 
average link G^, displayed in Figure 2a, the approximate method agrees well with 
the 11th order N = oo strong coupling series of ref. ||^ for (3 < 0.25. It also agrees 
the weak coupling N = oo series, 1 — 1/ (8/3) — 1/ (256/9^) + ... for large P, but in the 
intermediate coupling region it is below both curves and below the U{2) Monte Carlo 
data ||T^. Since the data for U{2) does not differ greatly from the data for SU{10) 



of |]14|, the dots in Figure 2a are probably close to the N = oo results. By allowing 
a and 6 to be a polynomial in Ga of sufficient high order, the approximate large N 
method can be adjusted to produce more terms in the strong coupling series. When 
this is done, agreement for the average link is obtained to less than 5% throughout 
the entire coupling constant region. 

From the mass gap, mA^^ = f{g), one can obtain a lattice beta function j3{g) 
(not to be confused with the inverse coupling) via —(3{g)/g = fl{g^^- In Figure 
2a, the approximate method is compared to the strong coupling result and to 
perturbation theory for which —f3{g)/g = g'^N/{16TT) + . . .. Excellent agreement is 
obtained for (3 < 0.25. In weak coupling, the approximate large N method gives 
results for —P{g)/g which are twice as large because the mass gap is exponentially 
suppressed only by half as much as the renormalization group result. 

VI. Conclusion 

In this work, we have used a vector-model-like method to linearized the lattice 
Schwinger-Dyson two-point function in the U (N) chiral model. For d <2, the approx- 
imate method produces no phase transition. This is in agreement with expectations 
for the finite case: The system is believed to be in the disordered phase and to 
have a mass gap for all couplings. For = oo, a phase transition occurs R P, p!6 



even though the system is disordered on both sides of the transition. For d < 2, 
our approximate method does not see this large transition, probably because the 
transition is rather mild and because it occurs in the intermediate coupling constant 
region, where the method is least accurate. For > 2, a phase transition occurs for 
all and is expected to be first order [0, |1^ . The approximate method correctly sees 
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the transition and correctly predicts the nature of the phases: symmetry breaking 
with massless Goldstone bosons at weaking couphng, while a disordered phase with a 
mass gap in strong coupling. Thus the method gives what-is-believed-to-be the cor- 
rect qualitative phase diagram in all dimensions for the finite case. By adjusting 
the method to agree with strong coupling expansions, the approximate method gives 
reasonably good results throughout the entire coupling constant regime for average 
link and mass gap. 

It is of interest to adapt the method to lattice gauge theories. The Schwinger- 



Dyson equations lead to relations among Wilson loops [T^, |T9|, |2^. A non-self- 
intersecting loop is related to loops modified at a link by the addition of plaquettes 
in the 2d directions minus loops modified by the addition of 2d "twisted plaquettes." 
The structure of the equation is similar to eq.(^ in the sense that the inverse cou- 
pling (3 multiplies the two terms and there is a sizeable cancellation between the two 
terms in weak coupling. The analog approximate large N method involves replacing 
each twisted-plaquette term by an untwisted loop and an original loop multiplied 
respectively by coefficients a and h. We have three comments on this approximate 
approach to large N for lattice gauge theories, (i) Although Wilson loops are related 
to Wilson loops after the substitution of the twisted-plaquette terms is made, there 
appears to be no way to find an analytic expression for Wilson loops when d ^ 2. 
This problem is probably related to the lack of an analytic description of free lattice 
string theory, (ii) When d > 2, the method appears to depend on the link on which 
the Schwinger-Dyson equations are derived. This implies that certain consistency 
issues must be addressed, (iii) When d = 2, the approximation can be consistently 
applied to the lattice gauge theory. For a U (N) gauge group, the results are identical 
to those of the U (N) chiral model. The value of a non-self-intersecting Wilson loop 
WL of area A is 

where A is given in eq. (|T3). The choice a = - {{TrUp)/NY and b = {{TrUp)/N) + 
{{TrUp)/Ny, where {TrUp)/N is the average plaquette (these values of a and b 
are the gauge-theory analog of eq. (|2^) ), leads to good results. The self-consistency 
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condition becomes {{TrUp)/Ny + {TrUp)/{N(3) = 1. The value of the Wilson loop 
then becomes 

{{TrUp)/Nf . (26) 

In the exact large N limit, the Wilson loop value is given by eg. (|26|) with (TrUp) /N = 
P for P < 0.5 and {TrUp)/N = 1 — 1/(4/5) for (3 > 0.5. Thus the discrepancy between 
the approximate and exact large N results is the same as in the d = 1 chiral model: 
The error is at most 6.7% near /3 = 0.6 and there is good agreement in the strong 
and weak coupling regimes. When a continuum limit is taken, the exact continuum 
solution is obtained. 

The picture of mass-gap generation for spin models in c? < 2 has an analogy for 
gauge models in d < 4. In the lattice gauge theory, the strong coupling expansion 
involves sums over closed surfaces. From a Euclidean viewpoint, such surfaces can 
be thought of as the propagation of strings. When two surfaces overlap, there are 
repulsive forces. In fact, the surface-surface interactions involve potentials which are 
exactly the same as in the chiral- model case 0. The quantum vacuum involves fluc- 
tuations in which strings appear, expand, contract and self-annihilate. This creates 
a medium in which any individual string must propagate. When d < 4, two large 
surfaces generically intersect, and since forces are repulsive, the medium creates a 
pressure which inhibits strings from becoming large. If this mechanism is sufficient 
robust as (3 increases and becomes smaller, then the string tension will remain 
non-zero in the continuum limit. The same mechanism can inhibit the spanning 
surface of a Wilson-loop computation from becoming large. Hence one arrives at a 
possible physical picture for confinement and mass-gap generation in c? < 4 gauge the- 
ories. The idea fails for d > 4 because two two-dimensional surfaces do not generally 
intersect. 
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Figure Captions 

Figure 1. The Average Link Ga Versus /3 for the N — oo d — 1 Case. 

Figure 2a. The Average Link Versus /3 for the N = (yo d = 2 Case. The dots 
are N = 2 Monte Carlo data. 

Figure 2b. The Beta Function —P{g)/g Versus the Inverse Couphng /3. The 
N — oo perturbative term comes from renormahzation group analysis. 
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